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ABSTRACT 

Although the giant radiogalaxy M 87 harbors many distinct regions of broad-band nonthermal 
emission, the recently reported fast variability of TeV 7-rays from M 87 on a timescale of days 
strongly constrains the range of speculations concerning the possible sites and scenarios of particle 
acceleration responsible for the observed TeV emission. A natural production site of this radiation 
is the immediate vicinity of the central supermassive mass black hole (BH). Because of the low 
bolometric luminosity, the nucleus of M 87 can be effectively transparent for 7-rays up to energy 
of 10 TeV, which makes this source an ideal laboratory for study of particle acceleration processes 
close to the BH event horizon. We critically analyse different possible radiation mechanisms in this 
region, and argue that the observed very high-energy 7-ray emission can be explained by the inverse 
Compton emission of ultrarelativistic electron-positron pairs produced through the development of 
an electromagnetic cascade in the BH magnetosphere. We demonstrate, through detailed numerical 
calculations of acceleration and radiation of electrons in the magnctospheric vacuum gap, that this 
"pulsar magnetosphere like" scenario can satisfactorily explain the main properties of TeV gamma-ray 
emission of M 87. 

Subject headings: gamma rays: theory — black hole physics — galaxies: active — galaxies: individual 
(M 87) 



1. INTRODUCTION 

M 87, a nearby giant radio galaxy, located at a dis- 
tance of d ~ 16 Mpc (Tonry 1991), hosts one of the 
most massive (M ~ 3 x 1O 9 M ) black holes (BH) in the 
nearby Universe (Marconi et al. 1997). M 87 contains a 
famous kpc-scale jet the high-resolution images of which 
detected at radio, optical and X-ray wavelengths show 
several prominent structures. Nonthermal processes play 
important, if not the dominant, role across the entire 
jet. The apparent synchrotron origin of the detected 
nonthermal emission, which extends from radio to X-ray 
bands, implies effective acceleration of electrons to multi- 
TeV energies. One may expect that protons and nuclei, 
which do not suffer radiative losses, are accelerated to 
much higher energies. Both the inner (sub-parsec) and 
large (kpc) scale parts of the jet of M 87 are possible 
sites of acceleration of protons to extremely high ener- 
gies (E ~ 10 20 eV). Therefore production of gamma-rays 
in different segments of the jet due to electromagnetic or 
hadronic processes is not only possible, but, in fact, un- 
avoidable. The jet of M 87 is observed at large angle, 
~ 20° (Biretta et al. 1999). Therefore, unlike blazars, we 
do not expect a strong Doppler boosting of the 7-ray flux. 
On the other hand, the nearby location of M 87 compen- 
sates this disadvantage (compared to blazars) and makes 
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several prominent knots and hot spots of the jet as po- 
tentially detectable TeV 7-ray emitters. 

In this regard, the discovery of a TeV 7-ray signal 
from M 87 by the HEGRA array of Cherenkov telescopes 
(Aharonian et al. 2003) and its confirmation by the HESS 
array of telescopes (Aharonian et al. 2006), was not a 
big surprise, especially given the rather modest apparent 
TeV 7-ray luminosity (few times 10 40 erg/cm 2 s). Sev- 
eral electronic and hadronic models have been suggested 
for explanation of TeV 7-ray emission of M87. The sug- 
gested sites of TeV 7-ray production range from large 
scale structures of the kpc jet (Stawarz et al. 2005) to 
a compact peculiar hot spot (the so-called HST-1 knot) 
at a distance 100 pc along the jet (Stawarz et al. 2006) 
and inner (sub-parsec) parts of the jet (Georganopoulos 
et al. 2005; Rcimer et al. 2005). 

While gamma-ray observations cannot provide images 
with an adequate resolution which would allow localisa- 
tion of sites of 7-ray production, the variability studies 
can discard or effectively constrain the suggested models. 

The continuous monitoring of M 87 with the HESS 
telescope array during the period 2003-2006 not only re- 
vealed statistically significant fluctuations of the TeV flux 
on a yearly basis, but, more excitingly, an evidence of 
fast variability on timescales of At ~ 2 days was found 
in the 2005 dataset (Aharonian et al. 2006). This re- 
quires a very compact region with a characteristic linear 
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size R < AtSj ~ 5 x 10 15 5j cm, where Sj is the Dopplcr 
factor of the relativistically moving source (throughout 
the paper we will use the system of units in which the 
speed of light c = 1). Note that since the mass of the 
BH in M 87 is well established (Marconi et al. 1997), 
the Schwarzschild radius is estimated quite accurately 
i? Schw = 2GM ~ 10 15 [M/3 x 10 9 M Q ] cm (G is the 
gravitational constant). The expected minimal variabil- 
ity time scale (light crossing time of the black hole) for 
a non-rotating BH is T\ s = 2i?g c h w / — 6 x 10 4 s~ 1 day, 
while it is two times less for the maximally rotating Kerr 
BH. The observed variability time scale of TeV emis- 
sion of ~ 2 days indicates that the emission originates 
within the last stable orbit of rotation around the black 
hole, unless the radiation is produced in relativistically 
moving outflow with a Doppler factor > 10. Although 
there are sound arguments against M 87 being a blazar 
(Fabian 2006), one cannot in principle exclude that at 
the base of its formation (close to the BH), where 7-rays 
are produced, the jet is pointed to us, and only later, it 
deviates from our line of sight. In this regard, both lep- 
tonic (Georganopoulos ct al. 2005) and hadronic (Rcimcr 
et al. 2005) models suggested for TeV radiation of M 87 
cannot be discarded. Note however, that both models 
predict rather steep energy spectra contrary to the ob- 
served hard 7-ray spectrum extending to E > 10 TeV 
(Aharonian et al. 2006). Whether these models are flex- 
ible enough to reproduce the observed TeV 7-ray spec- 
trum by tuning the relevant model parameters and in- 
troducing additional assumptions, should show further 
detailed theoretical studies. 

Finally, one should mention that formally one may as- 
sume that 7-rays are produced in a compact region far 
from the central engine. In this regard, the famous HST- 
1 knot has certain attractive features which make this 
compact structure a potential site of particle acceleration 
and 7-ray production (Stawarzet al. 2006). Although the 
favored size of this structure is in the range between 0.1 
and 1 pc, and thus contradicts the observed TeV 7-ray 
variability the lack of robust lower limits on the size of 
HST-1 leaves the 7-ray production in this peculiar knot 
as a possible option. Nevertheless one should note that 
the location of HST-1 at a distance of 100 pc from the 
central BH requires (almost) unrealistically tight colli- 
mation of the jet. 

In this paper we assume that the TeV 7-ray produc- 
tion takes place close to the event horizon of the central 
supcrmassivc BH, and show that the acceleration of elec- 
trons in a vacuum gap in BH magnctosphere can explain 
the general characteristics of the TeV 7-ray emission ob- 
served from M87. 

Such mechanism of 7-ray emission from the vicinity of 
a black hole is a close analog of the mechanism of pulsed 
7-ray emission from the vicinity of neutron stars in pul- 
sars. The similarity between the electrodynamics of the 
pulsar and black hole magnetospheres was discussed in 
the seminal paper of Blandford & Znajek (1977) in the 
context of pair production and energy transfer of rota- 
tional energy of a black hole through the Poynting flux. 
The mechanisms of emission of high energy 7-rays from 
the direct vicinity of black holes by electrons and protons 
accelerated in the electric field of vacuum gaps were dis- 
cussed by Blandford & Znajek (1977) (electron curvature 
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Fig. 1. — Top panel: optical depth for 7-rays produced in the 
vicinity of the black hole in the cases of the infrared source of size 
50-Rschw (black solid curve) and i?Schw/2 (grey dashed curve). The 
estimate of (Cheung et al. 2007) is shown by a cross. The spec- 
trum of the infrared background in the source is the one shown 
by red dashed curve in Fig. 8). Bottom panel: attenuation of 
7-rays in M 87 due to photon-photon pair production. The in- 
ternal absorption (thin solid curve) is dominated by interactions 
with the infrared radiation of the compact source in the core of 
M 87. The external absorption due to the interaction with the 
diffuse radiation fields within the elliptical galaxy M 87, the 2.7 
K CMBR and the diffuse extragalactic infrared background pho- 
tons leads to a further suppression of the 7-ray flux, shown by thick 
solid curve. Attenuation in the case when the 7-ray emission is dis- 
tributed throughout the infrared source is shown by thin (intrinsic 
absorption in the source) and thick (modification during the prop- 
agation through the galactic and extragalactic background light) 
dashed curves. Dashed grey curve shows the absorption of the 7- 
ray flux in the case of a "maximally compact" infrared source of 
the size Rm = -Rschw/2 (tee text). 

radiation), Bcskm et al. (1992) (inverse Compton scatter- 
ing by electrons) and Levinson (2000) (proton curvature 
radiation) . 

2. INTERNAL ABSORPTION OF 7-RAYS 

The observed infrared luminosity of the nucleus of 
M 87, vL v ~ 10 40 - 41 erg/s (Perlman et al. 2001; 
Whysong & Antonucci 2004) is 6.5 orders of magnitude 
lower than the Eddington luminosity of a 3 x 10 9 Af Q BH. 
In this regard the BH of in M 87 is similar to the super- 
massive BH in the center of our galaxy. In both cases 
the low bolometric luminosity of the nucleus makes the 
"central engine" of activity, i.e. the vicinity of the event 
horizon of the supermassive BH, transparent to the very 
high energy (VHE) 7-rays (Aharonian & Ncronov 2005). 

In the isotropic field of background photons, the cross- 
section of photon-photon pair production depends on the 
product of energies of colliding photons, s — Ee/m^- 
Starting from the threshold at s = 1, the cross-section 
ct 77 rapidly increases achieving the maximum <tq ~ 
<7t/5 ~ 1.3 x 10 -25 cm 2 at s ps 4, and then decreases 
as s^Hns. Because of relatively narrow distribution of 
<7 77 (s), gamma-rays interact most effectively with the 
infrared background photons of energy 

e ps l(E/l TeV) -1 eV . (1) 

Thus the optical depth for a gamma-ray of energy E 
produced in the center of the infrared source of the size 
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Rm and the luminosity Lm at energy given by Eq.(l) 
can be written in the form 



t(E,R ir ) 



LlRCTr, 
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The dependence of the optical depth on gamma-ray en- 
ergy is determined by the spectral form of background ra- 
diation n(e) = Lm(e) / (At: Re) . In particular, in the case 
of power-law spectrum with photon index T (n(e) oc e~ r , 
or L m (e) oc e~ r+1 , one has t(E, Rm) oc E v ~ x . Ac- 
curate numerical calculations of the optical depth for 
the spectral energy distribution of the compact infrared 
source in the nucleus of M87 (see Fig. 8, Section 5) and 
normalized to the source size Rm = 50i?schw is shown 
by black solid curve in the upper panel of Fig.l. Since 
t(E,Rir) oc -Rjr 1 , the optical depth does not exceed 1 
even at the highest detected energies of gamma-rays of 
about 10 TeV, provided that infrared source is larger 
than 50i?schw This is demonstrated in the lower panel 
of Fig. 1. Note that the recent claim by Cheung et al. 
(2007) that the central region of M87 is excluded as a 
site of the TeV emission because of absorption of 7-rays, 
is misleading. The authors obtained very large optical 
depth relevant to the energy ~ 25 TeV and assuming an 
extremely compact infrared source with a linear size of 
R = R g = i?schw/2 (the estimate of the optical depth 
by Cheung et al. (2007) is shown by a cross in the upper 
panel of Fig. 1, and the dependence of the optical depth 
on the 7-ray energy is shown by the grey dashed curve). 
Although formally one cannot rule out such a compact 
size of the IR source, a significantly larger size cannot be 
a priory excluded either. Moreover, there are not special 
reasons to assume that the infrared source is located very 
close to the event horizon. 

For the nucleus of M 87, there are no direct mea- 
surements of the size of the infrared source. Observa- 
tions in the microwave band at 43 GHz suggest that the 
size of the source at the mm wavelength is limited by 
5 x 10 16 cm, or approximately 50i?schw Lower angu- 
lar resolution of the infrared telescopes does not allow 
us to constrain (or marginally resolve, see Perlman et 
al. (2001)) the size of the nuclear source to < 10 pc 
(Whysong & Antonucci 2004). However, even assuming 
that the size of the infrared source is comparable to the 
size of the microwave source, one finds from the above 
estimate that the nucleus can be transparent to 7-rays 
with energies up to ~ 10 TeV. 

Even in the case of the "maximally compact" infrared 
source of the size Rm ~ i?schw/2, the source is partially 
transparent for 7-rays. In spite of the fact that the opti- 
cal depth for 7-rays produced in the center of the infrared 
source 7-rays becomes very large at energies E > 10 TeV 
(see the dashed grey curve on the upper panel of Fig. 1), 
there is no catastrophic absorption of multi-TeV 7-rays. 
The reason is that in this case the source(s) of 7-ray 
photons are distributed throughout the infrared source 
and the thickness H of the transparent surface layer of 
the source is determined from the condition t(E, H) ~ 1. 
Assuming a homogeneous 7-ray source, one can find that 
the luminosity of the last transparent layer (r < 1) is only 
moderately, by a factor of H/Rm ~ t(E, Rm) (rather 



than by a factor of cxp(— t(E, Rm)) lower than the total 
luminosity of the source. The attenuation of the 7-ray 
flux in the case of the 7-ray source distributed through- 
out the infrared source is shown by dashed curves in the 
lower panel of Fig. 1 for both cases of 50i?schw (black) 
and i?Schw/2 (grey) size of the infrared source. We as- 
sume the suppression factor (1 + t(E, i?iR)) _1 which is 
an interpolation between no supression in the r = limit 
and 1/r suppression in the r> 1 limit. 

The 7-rays after they escape the nucleus are further at- 
tenuated due the pair production in the radiation fields 
both inside and outside the elliptical galaxy M87. The 
spectrum of emission from the galactic bulge of M 87 
sharply peaks at photon energies around ebuigc — 1 eV. 
Interactions of nuclear 7-rays with the photon back- 
ground in M 87 galaxy should therefore lead to maxi- 
mum absorption at 7-ray energies around E 1 ~ 1 TeV. 
The column density of infrared/optical photons in the 
bulge of the size i?buige and luminosity Lbuige along the 
line of sight is estimated as 
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which allows to estimate the optical depth of gamma-rays 
at 1 TeV: 



t m 87 (1 TeV) = a N ph {\ eV) 



(4) 
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Propagation of 7-rays through the cosmic microwave and 
infrared backgrounds over the way from M 87 to the ob- 
server leads to further absorption of the highest energy 
quanta. Photons with energies above 10 15 eV are com- 
pletely absorbed due to interactions with the 2.7 K mi- 
crowave background (the minimal propagation distance 
is ~ 8 kpc). Photons with energies above 100 TeV inter- 
act most efficiently with the far-infrared background pho- 
tons whose density is some 3 orders of magnitude lower 
than the density of the microwave photons. However, the 
mean free path of E > 10 TeV 7-rays interacting with 
far-infrared background, is still shorter than the distance 
to M 87 (16 Mpc). 

3. PHYSICAL PARAMETERS OF THE CENTRAL ENGINE 

The high resolution observations of the nucleus of M 87 
in X-rays with the Chandra observatory provide an im- 
portant information about the accretion onto the super- 
massive BH (Di Matteo et al. 2003). In particular, they 
give an estimate of the electron density of plasma with a 
temperature kT ~ 1 keV, 

n e ~ 0.1 cm" 3 (5) 

at the distance of the order of Bondi accretion radius, 
-RBondi — 5 x 10 5 i?schw The corresponding accretion 
rate inferred from this estimate is 



0.1M Q yr" 



(6) 



Interestingly, the observed bolomctric luminosity of the 
nucleus of M 87 is 4 orders of magnitude below the ex- 
pected nuclear luminosity corresponding to this accretion 
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0.1M Q /yr 



erg/s. (7) 



where r\ ~ 0.1 is the efficiency of conversion of the rest 
energy of accreting particles into radiation. This indi- 
cates that cither the accretion proceeds in a radiatively 
inefficient way, or the actual accretion rate is still lower 
than the one inferred from X-ray observations. 

To estimate the plasma density close to the event hori- 
zon of the black hole, one has to assume a certain radial 
density profile, n(r) ~ r -7 . Depending on the model of 
accretion flow, the index 7 can vary between 1/2 (this 
value is, in fact, a lower limit which can be realized for 
collisionlcss motions of individual particles in the cen- 
tral gravitational field) and 3/2. The lack of information 
about the accretion regime leads to a significant uncer- 
tainty of the plasma density near the event horizon, 
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(8) 



Regardless of the uncertainty of this estimate, one may 
conclude that the strength of magnetic field in the vicin- 
ity of the BH can not be very high. Indeed, assuming that 
the magnetic field is generated by the accreting matter, 
one can find that the energy density of magnetic field can 
not exceed the density of the total kinetic energy stored 
in the particles of the accretion flow. In this case even 
if the accreting matter moves with relativistic speed, the 
estimate of maximal possible magnetic field is (assuming 
that the matter density is n ~ n, max ) 



(87m,, 



; ) 1/2 ~ 10 G. 



(9) 



Thus, particle acceleration close to the BH horizon pro- 
ceeds in the relatively low-density and low-magnetic field 
environment which significantly limits the range of pos- 
sible mechanisms of VHE 7-ray emission. Even for the 
maximally possible acceleration rate, dE/dt ~ eB eq , one 
can find that particles accelerated in a region of a lin- 
ear size of about the Schwarzschild radius can not reach 
energies higher than 



-Emax < eB cq Rs 



chw 



10 18 eV, 



(10) 



unless the magnetic field is significantly larger than the 
cquipartition estimate, given by Eq. (9). 

4. GAMMA-RAY EMISSION FROM ACCELERATED 
PROTONS. 

Protons accelerated near the BH horizon can produce 
7-ray emission in the VHE band through several radi- 
ation mechanisms. For example, TeV emission can be 
synchrotron or curvature 7-ray emission which accom- 
panies proton acceleration (Levinson 2000; Aharonian et 
al. 2002; Neronov et al. 2005). The energy loss time 
for protons emitting synchrotron radiation at the energy 
Esynch.p is short enough to explain the observed day-scale 
variability of the signal, 
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However, the energy of synchrotron and/or curvature 
photons produced by protons accelerated to the energy 



given by Eq. 
sion at 1-10 TeV, 



(10) is too low to explain the emis- 
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(12) 



(13) 



(assuming that typical curvature radius of proton tra- 
jectories is i?curv ~ -Rschw)- The 7-ray emission from 
the accelerated protons is thus expected in the 10 MeV - 
10 GeV energy region observable by GLAST, rather than 
in the TeV region visible by HESS. 

It is, in principle, not excluded that during short 
episodes of enhanced accretion the magnetic field can 
rise up to 10 3 G, which would, in principle, allow pro- 
ton acceleration up to the energies E max ~ 10 20 eV. 
Thus, the energy of curvature emission given by Eq.(13) 
can extend up to 10 TeV. Note, however, that even in 
this case the observed emission can not be related to 
the proton synchrotron radiation which has an intrin- 
sic self-regulated synchrotron cut-off at e sy nch < 0.3 TeV 
(Aharonian 2000), if the region of proton acceleration is 
spatially coincident with the region of synchrotron emis- 
sion. A potential problem of assumption about tran- 
sient enhancement (by 4 orders of magnitude, to produce 
the necessary increase of equipartition magnetic field, see 
Eq.(9)) of accretion rate is that it should result, in gen- 
eral, in a broad-band flaring activity of the nucleus of 
M 87, which, however, is not observed. 

VHE 7-rays are produced also in proton-proton (jyp) 
collisions. However, the interaction time of high-energy 
protons propagating through the low density medium 
(n < 10 7 cm~ 3 ) is quite large, 
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(14) 



~ 10~ 26 cm 2 is the proton-proton interaction cross- 
section) . Even if pp interactions would significantly con- 
tribute to the overall VHE 7-ray emission, they cannot 
explain the fast day-scale 7-ray flux variability of M 87. 
Therefore the observed variability should be referred to 
fast changes in the concentration of multi-TeV protons 
in the source, i.e. due to adiabatic or escape losses 
of protons, on timescales comparable to t vaT . The fast 
non-radiative losses versus slow rates of 7-ray produc- 
tion at pp collisions implies very low efficiency of conver- 
sion of the energy of parent protons to the VHE 7-rays, 
n = t-var/tpp < 3 x 10~ 4 . Thus to explain the 7-ray lu- 
minosity L 1 ~ 3 x 10 40 erg/s, the proton acceleration 
power should exceed k~ 1 L 1 ~ 10 44 erg/s which is just 
about the luminosity LBondi given by Eq. (7), assum- 
ing a conventional, 10% or so, efficiency of conversion of 
the rest mass energy of accreting particles into radiation. 
(Here we ignore a possible formation of gamma-rays in 
a relativistic outflow moving towards the observer which 
in principle would reduce by an order of magnitude this 
requirement). 

The energy requirements to the proton acceleration 
power can be somewhat relaxed if one invokes interac- 
tions of protons with the the surrounding radiation fields. 
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Although TeV 7-rays can be produced in a two-step 
process which includes Bethe-Heitler pair productions 
(P7 ~* P e+ i e ~) an d synchrotron radiation of secondary 
electrons, p"f interactions become efficient when they pro- 
ceed through the photomeson production channel. In or- 
der to interact with the photons of the infrared source 
with average energy eiR ~ 10 -2 eV, protons should be 
accelerated to E p ~ [200 MeV/em] m p ~ 2 x 10 19 eV. 
The number density of IR photons in the compact in- 
frared source is 
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For the average cross-section of the photo-pion produc- 
tion cross-section, a pl ~ 10~ 28 cm 2 , the interaction time 
of protons with infrared photons is 
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If the infrared source is very compact, _Z?ir ~ iischw/2, 
and the the accretion rate is transiently increased by 2 
orders of magnitude (to allow an order-of-magnitude in- 
crease in equipartition magnetic field and, as a conse- 
quence, proton acceleration to E > 10 19 cV), the p*y 
cooling time can be as short as the observed TeV vari- 
ability time scale. 

Note that the hypothesis of TeV 7-ray emission of M87 
based on the assumption of pj interactions, requires a 
very compact IR source with a size Rir ~ Rschw This 
implies strong absorption of gamma-rays with fast mul- 
tiplication of electron-positron pairs via Klein-Nishina 
cascades. Actually, photon-photon pair production is an 
important element of any p^ model; the observed spec- 
trum of TeV 7-rays cannot be explained by first gen- 
eration of ultra-high energy (> 10 15 ) photons from n° 
decays, and therefore requires production of secondary 
electrons which would provide broad-band emission in 
the TeV energy band. On the other hand, the copious 
pair production may lead to neutralization of the large 
scale (> i?schw) electric field, and thus to significant re- 
duction of the maximum achievable energy of protons 
given by Eq.(10). Since the rate of photomeson processes 
in the nucleus of M87 is very sensitive to the energy of 
protons, namely, it requires E p > 10 18 , the generation 
of large amount of secondary electrons may result in a 
dramatic drop of the rate of photo-meson production. 
A non-negligible contribution to the secondary electrons 
may come also from the Bethe-Heitler pair production, 
especially when the efficiency of photomeson production 
is suppressed (the energy threshold of this process is two 
orders of magnitude smaller than the energy threshold of 
the photomeson production). Whether this mechanism 
can explain the observed spectral and temporal charac- 
teristics of TeV 7-ray emission from M87, is a question 
which needs detailed numerical calculations. In any case 
it is clear that pj models can provide adequate efficiency 
only in the case of a very compact IR source with a size 
close to the Schwarzschild radius. 

5. GAMMA-RAY EMISSION FROM ACCELERATED 
ELECTRONS. 



5.1. Order-of-magnitude estimates 

The tough requirements of acceleration of protons to 
ultrahigh energies (E > 10 18 eV), as well as the rele- 
vant long cooling times challenge any interpretation of 
the day-scale variability of TeV 7-rays in terms of in- 
teractions of high-energy protons. The models based on 
acceleration of electrons do not face such problems, and 
arc likely to be responsible for the observed TeV 7-ray 
emission. 

The main emission mechanisms by electrons in 
the vicinity of the supermassive BH are syn- 
chrotron/curvature radiation and inverse Compton (IC) 
scattering. Electrons can be accelerated to multi-TeV 
energies only if the strength of the chaotic component 
of the magnetic field, B Tan ^, in the acceleration region is 
not too high. Assuming that electrons are accelerated at 
a rate dE/dt ~ neB m <i (k < 1 and -B rd is the ordered 
component of the magnetic field), from the balance of 
the acceleration and synchrotron energy loss rates one 
finds 
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Thus, even in the case of maximum possible acceleration 
rate (k = 1) electrons cannot emit in the 10-100 TeV 
band unless 

-Brand 5s 1 G . (18) 

In the ordered field the energy dissipation of electrons is 
reduced to curvature radiation loses. From the balance 
between the curvature loss rate and the acceleration rate, 
assuming that the typical curvature radius R curv of mag- 
netic field is comparable to the gravitational radius, one 
finds 
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Thus if the energy losses of electrons are dominated by 
curvature radiation, the maximum energy of accelerated 
electrons only weakly depends on the strength of the 
magnetic field. 

The IC loss rate is determined by the energy density of 
infrared radiation in the nucleus, U p h = Lir/ (47r i? 2 R c). 
The condition of the balance between IC loss rate and 
the electron acceleration rate gives 

W e Blf d R m ^ 



E e 
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(20) 
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k 1 / 2 eV . 



Note that this estimate is obtained assuming that the IC 
scattering takes place in the Thompson regime. However, 
highest energy electrons upscatter the infrared/optical 
radiation in the Klein-Nishina regime in which the effi- 
ciency of the IC scattering is reduced. A proper account 
of the decrease of the IC loss efficiency would result in 
higher electron energies exceeding the estimate of Eq. 
(20). 
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5.2. Electron acceleration in the vacuum gaps of BH 
magneto sphere 

So far we did not specify the particular mechanism of 
particle acceleration. In principle, several mechanisms 
can be responsible for the electron acceleration, but an 
obvious requirement which follows from the above esti- 
mates is that the "efficiency" parameter k for the accel- 
eration rate should be close to one, otherwise electrons 
would not reach the ~ 100 TcV energies, as it follows 
from Eqs. (17), (19), (20). Also, the irregular component 
of the magnetic field should not exceed 1 G. 

Large scale ordered electric fields, induced by rotation 
of black hole, are known to be responsible for particle ac- 
celeration and high-energy radiation in pulsars. A simi- 
lar mechanism of generation of large electric fields can be 
realised in the vicinity of a rotating BH placed in an ex- 
ternal magnetic field (Wald 1974; Bicak et al. 1976). In 
the case of pulsars, it is known that the force- free magne- 
tosphere possesses so-called "vacuum gaps" in which the 
rotation-induced electric field is not neutralized by redis- 
tribution of charges. The vacuum gaps work as powerful 
particle accelerators and sources of pulsed high-energy 7- 
ray emission. Vacuum gaps with strong rotation-induced 
electric field can be present also in the vicinity of a ro- 
tating black hole (Blandford & Znajck 1977; Beskin et 
al. 1992). Below we explore whether the observed VHE 
7-ray emission from M 87 can be explained by the emis- 
sion from the vacuum gaps formed in the magnctosphcrc 
of the supermassive black hole in M87. 

5.2.1. The magnetosphere of rotation-powered black hole 

Throughout the magnetosphere the component of elec- 
tric field directed along the magnetic field lines is neu- 
tralized by the charge redistribution, so that a force-free 
condition BLE is satisfied. The characteristic charge 
density needed to neutralize the parallel component of 
electric field in the magnetosphere of a BH rotating with 
an angular velocity f2 placed in an external magnetic field 
B is the so-called " Goldreich- Julian" density (Goldreich 
& Julian 1969) 



O ■ B Q1 
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aB n 
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(21) 



(a = il(GM) 2 , the BH rotation moment per unit mass, 
< a < GM, is a commonly used parameter of the Kerr 
metric describing the space-time of rotating black hole; 
see Appendix). 

In general, the charge distribution in the magneto- 
sphere is not static - additional free charges should be 
continuously supplied throughout the magnetosphere to 
compensate for the charge loss due to the magnctohydro- 
dynamical outflow. The inefficiency of charge supply can 
lead to the formation of " gaps" in the magnetosphere in 
which the parallel component of electric field is not zero 
and conditions for particle acceleration exist. 

In the case of pulsars, there are several potential ways 
to supply charged particles to the magnetosphere. First 
of all, the charge can be extracted directly from the sur- 
face of the neutron star. Electrons and positrons can be 
generated also due to pair production in very strong mag- 
netic field. Finally, electron-positron pairs can be created 
at interactions of 7-rays with low energy photons. 



Apart from the extraction of free charges from the sur- 
face of the compact object, the same mechanisms can 
in principle, be responsible for the charge supply to the 
magnctosphcrc in the case of black holes. However, in 
the particular case of the black hole in M87, the pair 
production of 7-rays in the magnetic field, (a mecha- 
nism, assumed e.g. in the Blandford & Znajek (1977) 
scenario) is not efficient because (1) the magnetic field 
cannot significantly exceed 10 G and (2) the energy of 7- 
rays emitted by accelerated particles cannot exceed 100 
TeV. On the other hand, the efficiency of the charge sup- 
ply via the pair production by 7-rays on the soft infrared 
background depends on the compactness of the infrared 
source (oc Lir/R). This process can be efficient only 
if 7-rays with energies above 10 TeV are present in the 
compact source. 

Since the 10 TeV 7-rays have to be produced by parti- 
cles accelerated to energies above 10 TeV, the gap(s), 
in which electric field component along the magnetic 
field lines is not neutralized by the charge redistribu- 
tion, should be present in the magnetosphere. In a self- 
consistent scenario the height of the gap(s) is limited by 
the condition that 7-rays emitted by the accelerated par- 
ticles do not produce e + e~ pairs within the gap. 

In order to estimate whether particle acceleration and 
high-energy emission from the vacuum gaps can be re- 
sponsible for the observed VHE luminosity of M 87 one 
has to estimate the total acceleration power output in 
the gap. In spite of the fact that the potential drop in 
the gap can be enough to accelerate charge particles to 
energies as high as 10 18 eV, strong radiative losses limit 
the maximum energy of electrons to 10 3 TcV. This means 
that the propagation of electrons through the gap pro- 
ceeds in a "loss-saturated" regime: all the work done 
by the gap's electric field is dissipated through the syn- 
chrotron/curvature and/or IC radiation. The rotation 
induced electric field near the BH horizon has a strength 
(see Appendix) £ ~ [a/GM]B OI( ±. For each electrons 
propagating in the gap, the energy loss rate is estimated 
as dE/dt ~ eS ~ e [a/{GM)] B orA } The density of elec- 
trons in the gap is limited by the Goldreich- Julian den- 
sity, given by Eq. (21). If the size of the infrared source 
is large enough, so that the gap height is not limited by 
pair production, the size of the gap is estimated to be 
about H ~ i?schw Taking into account that the volume 
of the gap is roughly R^ c h w H, the total number of elec- 
trons in the gap can be estimated as R'$ C T Ilvl Hn q . Then 
the total power output of the gap is 
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(22) 
erg/s 



Thus, if the angular momentum of the black hole is 
large enough, the nonthermal power of the vacuum gap 
can be as large as the observed TcV gamma-ray luminos- 
ity. 

5.2.2. Numerical modelling of acceleration and radiation of 
electrons in the gap of the black hole magnetosphere 



1 This implies that the acceleration efficiency re is re ~ [a/ (GM)}. 
For the extreme rotating black hole with a = GM, the acceleration 
reaches the maximum possible rate with re ~ 1. 
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Location of the gaps in the BH magnctosphcre de- 
pends on the structure of both the accretion flow and 
the magnetic field near the event horizon. In this re- 
gard, it should be noted that even after four decades of 
intensive theoretical study of physics of pulsar magne- 
tospheres, the details of the geometry of vacuum gaps 
remain uncertain. Nevertheless, the basic properties of 
particle accelerators operating in the vacuum gaps can be 
understood with a reasonable accuracy and confidence. 

We have developed a numerical code which allows, for 
the given geometry of the gap and configuration of the 
magnetic field, a quantitative study of energy distribu- 
tions of electrons accelerated in the vacuum gap and as- 
sociated electromagnetic radiation. For demonstration 
of the importance of the VHE 7-ray emission from the 
vacuum gaps, in this paper we have chosen a simple ge- 
ometry of the gap, namely we assumed that the gap oc- 
cupies a spherical layer above the BH horizon and has the 
height of about the size of the event horizon. The geome- 
try of electromagnetic field is assumed to be given by the 
solution of Maxwellian equations in Kerr metric, which 
corresponds to an asymptotically constant magnetic field 
inclined at an angle \ with respect to the rotation axis of 
the black hole (Bicak et al. 1976). The analytical solu- 
tion of Maxwellian equations are given in Appendix. The 
initial locations of electrons are assumed to be homoge- 
neously distributed either throughout the gap or over the 
boundary of the gap. The initial energies of electrons are 
assumed to be equal to the rest energy. Trajectories have 
been numerically integrated taking into account effects 
of General Relativity and energy losses of electrons due 
to synchrotron/curvature radiation and inverse Comp- 
ton scattering (see Appendix for technical details). The 
spectra and angular distributions of the synchrotron and 
curvature radiation are calculated by tracing the photon 
trajectories through the Kerr space-time metric from the 
emission point to infinity. 

As an example of numerical modelling, in Fig. 2 
we show distributions of average energies of electrons 
propagating in the spherical layer occupied by the gap. 
The magnetic field is assumed to be inclined at an an- 
gle x — 20° with respect to the rotation axis. The 
two left panels show angular distributions of the aver- 
age energy of individual electrons and the power of syn- 
chrotron/curvature emission, as measured in the Zero 
Angular Momentum Observer (ZAMO) frame (Bardeen 
et al. 1972), close to the black hole horizon (see Appendix 
for details). One can see that the maximum energies of 
photons are achieved in two oppositely situated hot-spots 
determined by the direction of magnetic field. At the 
same time, the total power of radiation does not strongly 
depend on the latitude and longitude coordinates. 

The " dark strips" (one along the equatorial plane and 
two snake-like dark strips above and below the equato- 
rial plane) are clearly recognizable in the left panels of 
Fig. 2. The drop of energies is explained by the specific 
configuration of electromagnetic field in these regions. 
Namely, the dark strips surround the so-called "force- 
free" surfaces at which the rotation-induced electric field 
E is orthogonal to the magnetic field _B or d • 

The right panel of Fig. 2 shows the angular distribu- 
tion of average photon energies and emissivity after the 
photon tracing to infinity through the Kerr space-time 
metric. One can see that polar "hot spots" become more 



pronounced, mostly because the photons emitted from 
equatorial regions (from the ergosphere) have larger red- 
shifts and a significant fraction of these photons is just 
absorbed by the black hole. 

Fig. 3 shows evolution of the shape of the polar hot 
spots with an increase of the inclination angle of the mag- 
netic field. One can see that the shape of the hot spots 
becomes wider and irregular. Also, with an increase of 
the inclination of the magnetic field the average photon 
energy decreases, which is explained by the fact that the 
acceleration of electrons is most efficient when the elec- 
tric field is aligned with the magnetic field. With the 
increase of the magnetic field inclination angle %i the 
regions of aligned electric and magnetic field (situated 
close to the rotation axis of the BH in the case \ = 0°) 
disappear. 

In Fig. 4 a typical spectrum of electrons accelerated 
in the spherical vacuum gap close to the BH horizon is 
shown. It is assumed that an extreme rotating BH (a = 
GM) is placed in 1 G magnetic field inclined at an angle 
of 60° with respect to the rotation axis. The size of 
the infrared emission region is assumed i?iR = 10iis c hw 
The three energy spectra shown in Fig. 4 correspond to 
different strengths of the random component of magnetic 
field. If the random magnetic field is smaller than 10 -3 
fraction of the ordered field, electrons propagating in the 
gap reach energies up to ~ 10 16 eV. With an increase of 
the random component of magnetic field, the synchrotron 
losses start to dominate which leads to reduction of the 
maximum electron energies energies, in a good agreement 
with the qualitative estimates of Section 5.1. 

5.3. Direct synchrotron/curvature and IC radiation 
from the acceleration process 

The radiative losses through both syn- 
chrotron/curvature and IC channels are released in 
the form of high energy 7-rays. The energy of Compton 
upscattered photons (in Thompson regime) is 



eic = 0.4 
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(23) 



The IC scattering of E e <; 10 TeV electrons on IR pho- 
tons proceeds in the Klein-Nishina regime, thus eic — E e . 
The curvature radiation peaks at significantly lower en- 
ergies, 
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(24) 

Since electron acceleration in the gap proceeds in the 
"loss saturated" regime, the calculations of the spectral 
and angular distributions of radiation accompaning the 
acceleration process, requires "self-consistent" approach 
in which the spectrum of radiation is calculated simulta- 
neously with the spectrum of parent electrons. The al- 
gorithm of self-consistent calculations used in this work 
is briefly described in the Appendix. 

Some results of self-consistent calculations of the 7-ray 
production spectra as functions of the viewing angle and 
the inclination of the magnetic field are shown in Figs. 5 
and 6, respectively. Fig. 5 demonstrates the difference of 
production spectra of 7-rays emitted along the direction 
of magnetic field (the region marked "on" in Fig. 2) and 
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Fig. 2. — Angular distributions of synchrotron/curvature photon energies (top) and 7-rayproduction rate (bottom) for emission from a 
spherical vacuum gap close to an extreme rotating BH placed in an magnetic field B = 1 G inclined at an angle \ = 20° with respect to 
the BH rotation axis. Two left panels show the angular distributions as seen in ZAMO frame close to the black hole horizon. Two right 
panels show the angular distributions after photon tracing to infinity. The scales are logarithmic and cover two decades from maximum 
(yellow on the top panels, red on the bottom panels) scale, max, down to the O.Olmax (black). Red and blue boxes on the bottom right 
panel show the regions used to extract spectra seen from "on hot spot" and "off hot spot" directions as shown in Fig. 5. 



away from this direction (the region marked "off' in Fig. 
2). Fig. 6 demonstrates the dependence of the 7-ray 
production spectra on the inclination angle of magnetic 
field. In both figures the low-energy (MeV-GeV) peak 
is due to the synchrotron/curvature radiation, while the 
high energy (TeV-PeV) peak is formed due to the IC 
scattering in the Klein-Nishina regime. 

5.4. Isotropic TeV emission from secondary 
pair-produced electrons 

The spectral energy distributions shown in Figs. 5 and 
6 correspond to the production rates of the first gener- 
ation 7-rays. They have essentially anisotropic distri- 
bution, thus the calculations of fluxes detected by an 
observer contain large uncertainties, mainly because of 
the poor knowledge of the source geometry. However, 
due to the internal and external absorption of > 10 
TeV 7-rays, the observer detects only a tiny fraction 
of the first generation 7-rays. While interactions with 
external photon fields lead to real attenuation of the 7- 
rayflux, the internal absorption is essentially recovered 
due to radiation of the pair produced electrons of sec- 
ond and further generations. Interestingly, the devel- 
opment of an electromagnetic cascade in radiation field 
of the infrared source may lead to "isotropisation" of 
the 7-ray source. Indeed, the absorption of first gener- 
ation 7-rays leads to deposition of e + e~ pairs through- 
out the infrared source volume, J£ir. If the latter is sig- 
nificantly larger than the volume corresponding to the 
vacuum gap (i.e. i?iR, 3> i?Schw)> the magnetic field in 
the IR source can be dominated by the irregular com- 
ponent which would effectively isotropisc the directions 



of secondary electrons. Correspondingly, the secondary 
radiation from e + e _ pairs will be emitted isotropically. 
For any reasonable magnetic field, the synchrotron radia- 
tion of secondary electrons is produced at energies signif- 
icantly below 1 TeV. Therefore for explanation of the ob- 
served TeV gamma-radiation one should assume that the 
energy losses of electrons are dominated by IC scatter- 
ing, i.e. the magnetic field in the infrared source should 
be significantly less than B = (L m /2Rf R ) 1 / 2 ~ 0.1 G. 
If so, the absorption of first generation gamma-rays will 
trigger an electromagnetic (Klein-Nishina) cascade. 

In Fig. 7 we show the resulting spectrum of gamma- 
radiation expected from the internal absorption of first 
generation gamma-rays. It consists of the isotropic com- 
ponent associated with the cascade in the infrared source 
(green dashed curve) and the primary anisotropic com- 
ponent whose intensity is uncertain since it strongly de- 
pends on the orientation of the observer with respect to 
the magnetic field direction. In Fig. 7 the thin solid 
red line corresponds to the sum of these two compo- 
nents, while the thick red solid line shows the result 
of absorption of the summary spectrum in the infrared 
source, in the elliptical galaxy M87 and in the intergalac- 
tic medium (see thick solid curve in Fig. 1). The curve 
is normalized to the observed flux of 7-rays at 0.5 TeV. 
The comparison of the calculated 7-ray spectrum shows 
quite a good agreement with the HESS measurements 
up to E ~ 10 — 20 TeV. One should note, however, that 
the agreement with the observations should not be over- 
emphasized, since we consider a "toy model" aimed to 
demonstrate the importance of TeV emission from the 
vacuum gaps in the magnetosphere. 
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Fig. 3. — Evolution of the shape of the polar hot spots with 
an increase of the inclination angle of magnetic field. The figures 
show the angular distribution of the energies of photons of syn- 
chrotron/curvature radiation traced to infinity. The color scale 
and parameters of numerical simulations are the same as in the 
top right panel of Fig. 2: maximum yellow corresponds to pho- 
ton energies 10 GeV, minimum (black) to photon energies below 
0.1 GeV. 



Finally in Fig. 8 we show the broad-band spectral en- 
ergy distribution (SED) of the resulting radiation and 
compare the model curve with observed fluxes of the nu- 
cleus of M87 at infrared, X-ray and TeV gamma-rays. 
Two broad peaks in the SED correspond to synchrotron 
radiation and inverse Compton scattering of secondary 
(cascade) electrons in the infrared source. The condition 
that the synchrotron emission from the infrared source 
should not exceed the observed flux in the X-ray band 
imposes an upper limit on the random magnetic field 
strength B < 0.1(i?/ij/i?schw) _1 G. The existing up- 
per limit on the M87 flux in the EGRET energy band 
(Sreekumar et al. 1994) imposes a restriction on the di- 
rect synchrotron/curvature emission from gap emitted in 
the direction of observer. 

Because of uncertainties of model parameters as well 
as the strong variability of radiation, we do not attempt 




Fig. 4. — Spectrum of electrons accelerated in the vacuum gap 
above the horizon of a maximally rotating BH of a mass M = 
3 X 10 9 Af Q placed in an ordered magnetic field B = 1 G inclined 
at 8 = 60° with respect to the BH rotation axis. Black solid line: 
the random magnetic field is 0.1% of the ordered one; blue dotted 
line: the random magnetic field is 1%; red dashed line: random 
magnetic field is 10%. 
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Fig. 5. — The production rate of first generation gamma-rays 
emitted by the spherical vacuum gap. The physical parameters 
are the same as in Fig. 2. Thick solid line: the total spectrum 
integrated over all directions. Dashed line: the spectrum integrated 
over the direction around the "hot spot" (the box marked "on" in 
Fig. 2). Thin solid line: the spectrum collected from the box 
marked "off" in Fig. 2. 
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Fig. 6. — The production rates of gamma-rays calculated for 
different values of the inclination angle of magnetic field. Solid 
line: \ = 20°; dashed line: \ = 60°; long-dashed line: \ = 90°. 
Physical parameters are same as in Fig. 2. 



to make a detailed fit to the broad band spectrum of 
the source, especially taking into account that the mea- 
surements at different energy bands correspond differ- 
ent epochs. On the other hand, the future simultane- 
ous studies of temporal and spectral properties of the 
broad-band emission can provide meaningful tests of the 
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Fig. 7. — Secondary emission from high-energy electrons in- 
jected into the compact infrared source in the nucleus of M 87 via 
the photon-photon pair production. Thick black solid line: om- 
nidirectional spectrum of primary emission from accelerated par- 
ticles; thin blue solid line: the primary spectrum only from the 
"off' direction (same as in Fig. 5); black and blue thin dashed 
lines show the spectra attenuated by the pair production in the in- 
frared source. Red dotted line shows the contribution of secondary 
cascade (isotropic) emission. 
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Fig. 8. — Spectral energy distribution of emission from the nu- 
cleus of M 87. Dashed blue line: isotropic synchrotron and IC emis- 
sion from secondary electron positron pairs injected via the pair 
production. Short-dashed green line: direct synchrotron/curvature 
and IC emission from electrons accelerated in the vacuum gap 
(strongly anisotropic, depends on the geometry of the vacuum gap). 
Solid blue line shows the total emission spectrum, which is the sum 
of the direct and cascade contributions. Dashed red line shows the 
model spectrum of soft photon background used for the calculation 
of IC scattering. This radiation component can come from a larger 
region and is not necessarily related to the particle acceleration in 
the vacuum gap. 



energy losses due to the curvature radiation and inverse 
Compton scattering require an extremely effective mech- 
anism of particle acceleration with a rate close to the 
maximum (theoretically possible) acceleration rate. 

In this paper we show that the observed TeV gamma- 
ray emission from M87 can be explained by electrons ac- 
celerated in strong rotation induced electric fields in the 
vacuum gaps in black hole magnetosphere. Generally, 
this model has many similarities with models of parti- 
cle acceleration in pulsar magnetosphercs. Our detailed 
modelling shows that the gamma-radiation from the cen- 
tral engine of M 87 consists of both first generation pho- 
tons emitted by particles accelerated in the gap (severely 
attenuated due to interactions with the internal and ex- 
ternal radiation fields) and second and further genera- 
tion (cascade) photons. If the first component dominates 
above 10 TeV, the cascade 7-rays contribute mainly to 
the < 10 TeV energy domain. The electron acceleration 
and 7-ray production in a very compact region close to 
the event horizon of the black hole naturally explains the 
observed variability of TeV 7-ray emission from M87. 
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proposed model and significantly reduce the relevant pa- 
rameter space. 

6. SUMMARY AND CONCLUSIONS. 

We have studied the mechanisms of production of vari- 
able TeV 7-ray emission from vicinity of the supcrmas- 
sive BH in the nucleus of M 87. Moderate accretion rate 
onto the black hole, inferred from the Chandra observa- 
tions, limits the magnetic field strength close to the black 
hole horizon - it cannot significantly exceed B < 10 G. 
This limits the maximum energy attainable by protons, 
E < 10 18 eV. None of the known mechanisms of 7-ray 
emission by protons can satisfactorily explain the ob- 
served temporal and spectral characteristics of the ob- 
served TeV 7-ray emission. On the other hand, severe 
radiative losses of electrons in a compact region close to 
the black hole significantly constrains the range of possi- 
ble acceleration scenarios. In particular the random com- 
ponent of the magnetic field cannot exceed 1 G. Thus, 
the acceleration takes place, most likely, in a region where 
the regular magnetic field significantly exceeds the ran- 
dom component of the field. Even so, the unavoidable 
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APPENDIX 

DETAILS OF THE NUMERICAL MODELLING OF PARTICLE ACCELERATION IN THE VACUUM GAPS IN 

MAGNETOSPHERES OF ROTATING BLACK HOLE. 

A self-consistent modelling of acceleration of electrons in the direct vicinity of event horizon of a BH requires (a) 
a full account of the effects of General Relativity and (b) a full account of the radiation reaction on particle motion, 
since electrons propagate most of the time in the "loss saturated" regime when the acceleration force is balanced by 
the radiation reaction force. Below we give some details of the modelling of trajectories electrons and photons in the 
vicinity of the black hole. 

The Kerr space-time. 

A rotating BH is described by two parameters: its mass M and the angular momentum per unit mass a < GM. 
The geometry of space-time in the vicinity of horizon is described by the Kerr metric 

ds 2 = -a 2 dt 2 + g lk [dx l + /3 l dt] [dx k + (3 k dt\ (Al) 

pVA p 2 2 £ 2 sin 2 

9rr = -r! gee = P ; 



^ — £ ' — A ' ' — p 2 ' 

2aGMr , 9 

P* = J A = r + a — 2GMr; 

E 2 = (r 2 +a 2 ) 2 - a 2 Asin 2 0; p 2 = r 2 + a 2 cos 2 (A2) 

The horizon is situated at r H = GM + y/(GM) 2 - a 2 . 

To understand the acceleration and energy losses of charged particles propagating close to the BH horizon, it is 
convenient to use an orthonormal (non-coordinate) frame 



£ d 2GMar d _ VA d Id __p d_ 



py/Adt Ypy/A d<t>' r p dr 1 9 pd6' * Esin0< 

carried by the so-called "zero angular momentum" observers (ZAMO) (Bardeen et al. 1972). The corresponding 

covariant basis vectors e l are given by 

pVA f p § i Esin0 2GMar sin . 

e u = dt; e = -^=dr; e" = pd6, e v = d(j> dt. (A4) 

£ VA p pE 



The electromagnetic field. 

In the reference frame (A3) the magnetic field inclined at angle \ with respect to the BH rotation axis is given by 
(Bicak et al. 1976) 

B f = r {Bu sin0cos0 [Ap 4 + 2GMr(r 4 - a 4 )} + 

£p 4 sm0 

B± [rcosi* - a sin ip] [p 4 (r sin 2 + GM cos 20)- 
GM(r 2 + a 2 )(r 2 cos 20 + a 2 cos 2 9)] } 

B § = - - j-j= {B||Asin6» [p 4 r + a 2 GM(r 2 - a 2 cos 2 9)(1 + cos 2 9)} + 

B±_ cos 9 [p 4 (( Ar - GMa 2 ) cos tp + a( A + GMr) sin i/)] - 
a 2 GMr sin 2 9 (r 2 (r - 2GM) + 2a 2 (r sin 2 9 + GM cos 2 9)) cosifj- 
a(A - 2GMr - 2a 2 cos 2 6) sin ib] } (A5) 



where 



tp = <t> H , In 

2 V / (GA/) 2 - a 2 



r - GM + y 7 {GM) 2 - a 2 



B«=B cosx, B ± = B smx- (A6) 



= l^fli A [2r 2 p 4 sin 2 9 - (£ 2 — 2GAfra 2 sin 2 0)(r 2 - a 2 cos 2 9)(1 + cos 2 0)] - 



r-GM- y/(GM) 2 

Rotation of the BH is responsible for the appearance of nonzero electric field whose components are 
aGM 
EAp 

B±r sin 9 cos 9 [2 [(VA - GMa 2 ) cosip - a(A + rGM) sin^] + 
(£ 2 - 2GMra 2 sin 2 0) [r 2 (r - 2GM) + 2a 2 (r sin 2 + GM cos 2 0) cos ^ 
o(A - 2GMr - 2a 2 cos 2 0) sin xb] ] } 

= aG l ! {2B|,rsin0cos0 [Ap 4 - (r 2 - a 2 ) [E 2 - 2GA/r(r 2 + a 2 )]] + 



B_L [2rp 4 [r(r sin 2 + GM cos 20) cos ip - a{r sin 2 + GM cos 2 0) sin %b] - 

(r 2 cos 20 + a 2 cos 2 0) [E 2 - 2GMr(r 2 + a 2 )] (a sin ip - r cos ib)] } (A7) 
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Equations of motion for a charged particle. 
The components of the four- velocity of a particle ?/' = dx^/dt in the orthonormal frame e- a (A4) are 

w' 1 = v h eldijdt (A8) 

where t is the time which would be locally by the ZAMO observers at a given point and t is the coordinate time which 
enters the metric (Al). For example, the <fr components of particle velocity in coordinate and orthonormal reference 
frame are related through 



£ sinf? 



d(j) 2GMar 



^dt £ 2 

The extra term f2 = IGMar /T, 2 is the angular velocity of the ZAMO frame at each point. From (A4) one can see that 



It is convenient to introduce a particle 7-factor in the orthonormal frame 7 = 1/y 1 — (v a )' 2 . 
Equations of motion in the orthonormal basis (A3) have the same form as in the flat space (Thorne et al. 1986) 

dz) -* -* " — * 

-£ = e(E + v x B) + m-yg + Hp + f rad (All) 
at 

where p is the particle momentum 

p & = mi v h (A12) 

g is the gravitational acceleration and H is the tensor of gravi-magnetic force. The force frad is the radiation reaction 
force. In the case of interest the time scales for acceleration in electromagnetic field and of the radiation reaction are 
orders of magnitude shorter than that of the motion of the particle in the gravitational field of the black hole. 

The radiation reaction force f ra d for the ultra-relativistic particles moving in external electromagnetic field is (see, 
e.g. (Landau & Lifshitz 1975)) 

frad = \(E + vx B) 2 -(v-(E + vx B)) 2 ] ^ (A13) 



3m 2 



Note that if particles move at large angle with respect to the magnetic field lines, this expression will describe mostly 
synchrotron energy loss. However, in the case when particles move almost along the magnetic field lines, the last 
equation will " mimic" the effect of curvature energy loss (taking into account the fact that the typical curvature radius 
of the magnetic field lines in the considered case is about i?schw)- 



REFERENCES 



Aharonian, F.A., 2000, New AR, 5, 377. 

Aharonian F.A. and Neronov, A. 2006, Ap.J., 619, 306 

Aharonian, F. A., Belyanin, A. A., Derishev, E. V., Kocharovsky, 

V. V., Kocharovsky, VI. V. 2002, Phys ReV. D, 66, id. 023005 
Aharonian, F. et al. (HEGRA collaboration) 2003, A&A, 403, LI 
Aharonian, F. et al. (HESS collaboration) 2006, Science, 314, 1424 
Bardeen, J.M., Press, W.H., Teukolsky, S.A., 1972, Ap.J. 178, 347. 
Beskin V.S., Istomin Ya.N., Par'ev V.I., 1992, Soviet Astronomy, 

36, 642. 

Bicak J., Dvorak L., 1976, Gen.Rel.Grav., 7, 959; see also Bicak J., 

Janis V., 1985, MNRAS, 212, 899. 
Biretta, J. A., Sparks, W. B., & Macchetto, F., 1999, ApJ, 520, 

621. 

Blandford R.D., Znajek R.L., 1977, MNRAS, 179, 433. 
Cheng K.S., Ho C, Ruderman M., 1986, Ap.J. 300, 500. 
Cheung C.C., Harris D.E., Stawarz L., 2007, Ap.J. Lett, accepted, 
arXiv:0705.2448. 

Di Matteo T., Allen S.W., Fabian A.C., Wilson A.S. Young A.J. 

2003, ApJ, 582, 133. 
Fabian, A. 2006, Science, 314, 1398 

Georganopoulos, M., Perlman, E. S., Kazanas, D. 2005, ApJ, 634, 
L33 

Goldreich P., Julian W.H., 1969, Ap.J. 157, 869. 

Harris D.E., Biretta A. J., Junor W., Perlman E.S., Sparks W.B., 

Wilson A.S., 2003, Ap.J., 586, L41. 
Heinz S.; Bcgelman M.C., 1997, Ap.J., 490, 653 



Landau L.D., Lifshitz, E.M., 1975 The Classical Theory of Fields, 

Oxford: Pcrgamon Press 
Levinson A., 2000, Phys. Rev. Lett., 85, 912. 

Marconi, A., Axon, D. J., Macchetto, F. D., Cappetti, A., Sparks, 

W. B., Crane, P. 1997, MNRAS, 289, L21 
Michel F.C., 2004, Ad.Sp.R., 33, 542. 

Neronov A., Semikoz D., Aharonian F., Kalashcv O., 2002, 

Phys.Rev. Lett., 89, 1101. 
Neronov A., Tinyakov P., Tkachev I., 2005, JETP, 100, 656. 
Perlman E.S., Sparks W.B., Radomski J., Packham C, Fisher R.S., 

Pina R., Biretta J.A. 2001, Ap.J., 561, L51. 
Reimer, A., Protheroe, R. J., Donea, A.-C. 2004, A&A, 419, 89 
Sreekumar P. et al., 1994, Ap.J. 426, 105. 

Stawarz, L., Siemiginowska, A., Ostrowski, M., Sikora, M. 2005 
Ap.J., 626, 120 

Stawarz, L., Aharonian, F., Kataoka, J., Ostrowski, M., 
Siemiginowska, A., Sikora, M. 2006, MNRAS, 370, 981 

Thorne K.S., Price R.H, Macdonald D.A., 1986, Black Holes: the 
Membrane Paradigm, Yale University Press. 

Tonry, J. L. 1991, ApJ, 373, LI 

Wald R.M., 1974, Phys.Rev., D10, 1680. 

Whysong D., Antonucci R., 2004, Ap.J., 602, 16. 

Young A.J., Wilson A.S., Mundell C.G., 2002, Ap.J., 579, 560. 



